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Method and associated device . for generating random 
numbers in a given range 

This disclosure is based upon French Application 
5 No. 0312435 filed October 24, 2003 and International 
Application No. PCT/FR2004/050510 , filed October 18, 
2004, the contents of which are incorporated herein by 
reference . 



10 BACKGROUND OF THE INVENTION 

The invention concerns a method of obtaining a 
random number between A and B from a generator 
producing random numbers lying between 0 and W-l, with 
N the size of the numbers produced by the generator, W- 

15 1 the maximum value taken by the random numbers 
produced, with for example W = 2 N and A, B any integer 
numbers, less than or greater than the number W. 

Such a situation occurs for example in an 
electronic component adapted to perform cryptographic 

2 0 calculations and comprising an N-bit random number 
generator, for example N = 8. The random numbers that 
it can produce are thus between 0 and W-l = 255, whilst 
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it would be desirable to have random numbers between 
for example 0 and 100 or between 300 and 10000. It 
should be noted that it suffices to determine numbers 
between 0 and 97 00 and then to add 3 00 to the number 
5 obtained in order finally to obtain a number between 
300 and 10000. 

Such a situation is found in practice in the 
majority of cryptographic applications, for example the 
DSA signature, the El Gamal signature or enciphering, 
10 the development of countermeasures against various 
attacks, etc. 

Several methods are already known for producing 
random numbers R between 0 and K from numbers between 0 
and W-l. These methods are in general implemented by 
15 software means used to control on the one hand a 
hardware generator that produces random numbers of size 
N and on the other hand calculation means performing in 
particular multiplication, addition, etc operations. 

A first known method comprises the following 

20 steps: 

a) determining the smallest integer number p such 
that K < WP - 1, 

b) producing p random numbers S 0 , S, . . . , S p _i and 

p-i 

forming the variable S = ^S i *W 

i=0 

25 c) if S > K, then returning to step b) , otherwise 

putting R = S 

R is the random number sought, between 0 and K. 

The equation S = /TS f *W l is a representation of the 

i=0 
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variable S decomposed/recomposed in base (W p_1 , . . . , W 1 , 
W°) . It would also be possible to note S 

Sp_iSp_ 2 . . . SiS 0 / a notation commonly used. 

A second known method comprises the following 

5 steps: 

a) determining the smallest integer number p such 
that K < WP - 1 , 

b) producing p random numbers S 0 , S, . .., S p _i and 

P-2 

forming the variable T^^S,. *PT and S = T + Sp-!*!/^" 1 

i=0 

10 c) if S > K, putting R = T otherwise putting R = 

S 

A third known method comprises the following 

steps : 

a) determining the smallest integer p such that K 
15 < WP - 1, 

b) producing p random numbers S 0 , S, S p -i and 

p-\ 

forming the variable S = ^ d S i *W i 

/=o 

c) putting R = S mod(K+l), that is to say the 
remainder of the whole-number division of S by K+l, 

20 also referred to as modular reduction of S by K+l. 

These three methods can be summarised by the 
following steps: 

a) producing p random numbers S 0 , S, . . . , S p _i, 
being the smallest integer number such that K <W P - 1 

2 5 and forming the variable S = Y j S i *W i 

b) determining the random number R from the 
variable S. 
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According to circumstances, during step b, R is 
obtained from S by repeating step b (first method) , 
taking account or not of the additional random number 
Sp_i (second method) or performing a modular reduction 
5 (third method) . 

It should be noted that, in the three methods, if 
a number between A and K+A is required, it suffices to 
add A to the number R obtained lying between 0 and K. 

The main drawback of the first method is a 

10 particularly long and especially unpredictable 
calculation time: the step of producing the p random 
numbers may be repeated numerous times without it being 
possible to predict at the start the number of 
repetitions of this step. 

15 The second and third methods have the main 

drawback of producing random numbers exhibiting a bias: 
amongst the numbers R produced in the range [0, K] , 
certain values are more probable than others. In other 
words, the numbers R produced are not perfectly random 

20 (non-uniform distribution) . This bias may have 

significant consequences on the security of the 
cryptographic systems liable to implement these methods. 
The security of cryptographic systems assumes in fact 
that the random numbers that they use are uniformly 

25 distributed (or at least close to a uniform distribution) 
in the range [0, K] or [A, K+A] wished for. 

Finally, the three methods are slow overall 
because they implement operations on large numbers, of 
size N (in the sense of the number of bits) greater than 

3 0 the size of the circuits used for the implementation. 
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This is because the number K in particular is any number 
and can be greater than W and therefore of size greater 
than N. The variable S can also be of large size. 
However, the implementation of operations on large 
5 numbers requires the implementation of complex methods 
expensive in terms of calculation time. 

DESCRIPTION OF THE INVENTION 

An essential object of the invention is to 
10 propose a method of constructing a random number R that 
is particularly rapid. 

Thus the invention proposes a cryptographic 
method during which use is made of a random number 
generator producing random numbers S± of size N fixed 
15 between 0 and W-l, with for example but not necessarily 
W = 2 N , in order to produce a random number R between 0 
and a predefined limiter K. 

The essential steps of a method according to the 
invention are as follows: 
20 E31: a random variable Si between 0 and W-l is 

produced, 

E32: if the random variable Si is strictly less 
than a coefficient Ki of the limiter K in base W, then 
the coefficient Ri of rank i of the random number R is 
25 equal to the random variable Si and then, for any rank J 
less than i, a random variable Sj between 0 and W-l is 
produced and Rj = Sj , 

E33: otherwise, if the said random variable is 
greater than the coefficient Ki of rank i of the limiter 
30 K in base W, then the said coefficient Ri is determined 
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from the random variable Si of rank i according to a 
predetermined function, and then the coefficient Ri_! is 
determined for the random number R of rank i-1 that is 
immediately lower by repeating steps E31 to E33 . 
5 Thus, in a method according to the invention, the 

coefficients Ri of the random number R required are 
sought one by one, commencing with the most significant 
coefficient R p -i. The physical generator of random 
numbers used thus produces random variables Si one by 

10 one, one variable at each iteration. 

In addition, the method is rapid since step E33 is 
executed a small number of times. This is because, as 
soon as one of the variables Si produced by the physical 
generator is less than the associated coefficient Ki of 

15 the limiter K, the method no longer requires the 
processing of the variables S j of rank less than i : thus 
a small number of coefficients of the number R, the most 
significant, are calculated the most often. 

Finally, compared with the known methods, a 

2 0 method according to the invention has the advantage of 
working on numbers of no more than N bits, N being the 
size of the registers and other calculation circuits of 
the devices used for implementation. For example, if W 
is equal to S N , the coefficients Ki resulting from the 

25 decomposition of K in base (W p " 1 , . . . W 1 , W°) are 
necessarily less than W and therefore with a size of no 
more than N bits. Likewise, the random variables Si 
produced by the physical random number generator are 
also of N bits. 
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to 



the 



essential 
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initialisation step and a step of recombination of the 
random number R, there are obtained: 



El: the limiter K is decomposed in base (W p 1 / W p 2 



. .., W°) (K = Y d K i *W i or K = K p " 2 ... K^K 0 ) , i being a 



loop index, Ki being a coefficient of the limiter K of 
rank i between 0 and W-l and p being the degree of the 
limiter K, 

E2 : a Boolean variable f is initialised to TRUE, 
E3 : the following operations are performed, in a 

loop indexed by i, i being an integer varying between 

p-1 and 0: 

E31: a random variable Si between 0 and W0-1 
is produced, 

E32 : if the random variable S± is strictly 
less than the coefficient Ki of rank i, then the 
Boolean variable f is set to FALSE, 

E33_l: if the random variable Si is strictly 
greater than the coefficient Ki of rank i and the 
Boolean variable f is TRUE, then the coefficient Ri 
of rank i is determined from the random variable Si 
of rank i according to a predefined function, 
E33_2: otherwise Ri = Si 
E34 : the loop indexed i is decremented, 
E4 : the random number R is determined by 
recombination of the random coefficients Ri in base 



W(i? = £i?,.*fT or R p -\ . . 




R X R°) . 



In concrete terms, as soon as the Boolean 
variable f is positioned at FALSE, it remains at this 
value since provision is not made for repositioning it 
at the value TRUE, except when E2 of the method is 
5 initialised. Step E32 is executed only if the variable 
f is TRUE; thus, as soon as the variable f is 
positioned at the value FALSE, step E3 3__l is no longer 
executed and the method according to the invention ends 
rapidly. 

10 A second objective of the invention is to propose 

a method of constructing random numbers whose 
distribution is uniform or can be made as close as 
desired to a uniform distribution. This objective is 
achieved by choosing a suitable function for the 

15 determination of the coefficient Ri from the random 
variable Si. 

According to a first embodiment of the method 
according to the invention, in order to determine the 
coefficient R± of rank i from the random variable Si of 
20 rank i (step E33_l) , the following substeps are 
performed : 

E33_ll: if the random variable Si is strictly 
greater than the coefficient Ki of the limiter K, then a 
new random variable Si is produced, 
25 E33_12 : step E33_ll is repeated until the random 

variable Si is less than the coefficient Ki of the 
limiter K, and then the coefficient Ri is equalised to 
the random variable Si. 

In such an embodiment, all the coefficients Ri 
3 0 obtained are numbers directly produced by the hardware 
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random number generator; and these coefficients are 
therefore perfect and the number R which results 
therefrom is also perfect. In other words the 

distribution obtained of the numbers R is uniform in 
5 the range [0, K] . 

According to a second embodiment, during step E33 
the coefficient Ri of rank i is chosen so as to be equal 
to part of the random variable Si, a part less than the 
coefficient K± . The said part corresponding in one 
10 example to a limited number of bits of the variable Si. 

According to a third embodiment, during step E33 
the random variable Si is reduced modulo Ki+1, the 
results of the reduction being the coefficient Ri 
sought . 

15 These latter two embodiments are rapid compared 

with the known methods, essentially because the work is 
done on small numbers. The distributions of random 
numbers obtained are however not uniform: the simple 
fact of truncating the variable Si or performing a 

20 reduction modulo Ki+1 necessarily introduces a bias. 

However, this bias is less compared with the methods of 
the prior art . 

Moreover, it is possible to reduce the bias of 
the methods according to the second and third 

25 embodiments proposed, as will be seen below. 

In a method according to the invention as 
described above, a random number R is constructed less 
than K from variables Si of size N produced by a 
perfectly random physical generator. The number R 
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obtained is biased, but the bias is small compared with 
a known method. 

For this, in the second or third embodiment, a 
coefficient Ri < Ki is constructed, in particular during 
step E33__l, from variables Si of size N. In order to 
reduce the bias introduced on the coefficient Ri, it is 
proposed to construct it using the same steps El to E3 
as for constructing the number R. In a sense, two 
similar methods are "interleaved" . This makes it 
possible to reduce further the size of the numbers on 
which the work is carried out, and consequently to 
reduce further the bias on the coefficient of R, and on 
the final number R. 

In concrete terms, in order to determine the 
coefficient Ri of rank i from the random variable Si of 
rank i (step E33_l) , steps El to E4 are executed using 
a base ((3 q_1 , P°) as the calculation base, P being 

an integer number strictly less than W and q being the 
degree of Ki in base p. 

Step E3 3 is thus broken down into the following 
substeps : 



E33_41: the coefficient Ki of rank i of the 



limiter K in base (P q "\ p°) ( K x = or K ± = 



(Ki) q _i . . . (Ki) i (Ki) o) / j being a loop index, (Ki)j being a 
number between 0 and P-l and q being a degree of the 
coefficient Ki, is decomposed, 

E33_42: a second Boolean variable g is 
initialised to TRUE, 
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E33_43 : the following operations are performed, 
in a loop indexed by j varying between q-1 and 0: 

E33_431: a random variable (Si)j between 0 
and P-l is produced, 
5 E33_432: if the random variable (Si)j is 

strictly less than the coefficient (Ki)j, then the 
second Boolean variable g is set to FALSE, 

E33_4331: if the random variable (Si)j is 
strictly greater than the coefficient (Ki)j and the 
10 second Boolean variable g is TRUE, then a 

coefficient (Ri) j is determined from the random 
variable (Si)j according to a predefined function, 
E33_ 4332: otherwise, (R±)j = (Sj.)j 
E33_434: the loop indexed j is decremented, 
15 E33_44 : the random number Ri is determined by 

recombination of the random coefficients (Ri)j in base (3 

(tf, or Ri = (Ri) q -i. . . (Ri)i(Ri)o) - 

As has just been seen above, by "interleaving" 
two methods, the bias of the random numbers R produced 

20 by the global method is reduced, whilst preserving a 
rapid global method. It is of course possible to 
imagine "interleaving" more than two methods, for 
example three or four, by decomposing, in step E3 3_4 3, 
the numbers in base y < p, and decomposing step E33_43 

25 in a succession of steps similar to steps E33_41 to 
E3 3__4 3 . 

In general terms, the more methods are 
"interleaved", the smaller the numbers on which the 
work is carried out: the duration of each step 
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decreases and the bias of the numbers produced by the 
global method also decreases. 

Another object of the invention is an electronic 
component adapted for implementing the method as 
described above. Such a component comprises in 

particular a generator producing random numbers of size 
N, and calculation circuits for performing operations 
on numbers of no more than N bits. 

According to the embodiment of the method to be 
implemented, the calculation circuits are adapted to 
perform operations of comparing two numbers, number 
truncation and modular reduction. 

The random number generator and the calculation 
circuits are preferably controlled by a software means 
stored in a memory of the component provided for this 
purpose . 

The invention also concerns a chip card 
comprising an electronic component as described above. 



